Quantum Reversibility and Echoes in Interacting Systems 
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In Echo experiments, imperfect time-reversal operations are performed on a subset of the total 
number of degrees of freedom. To capture the physics of these experiments, we introduce a partial 
fidelity Mb((), the Boltzmann echo, where only part of the system's degrees of freedom can be 
time-reversed. We present a semiclassical calculation of Mb(£). We show that, as the time-reversal 
operation is performed more and more accurately, the decay rate of Mb (t) saturates at a value given 
by the decoherence rate of the controlled degrees of freedom due to their coupling to uncontrolled 
ones. We connect these results with NMR spin echo experiments. 



PACS numbers: 05.45.Mt,03.65.Ud,05.70.Ln,03.67.-a 

One of the central problems faced by the founders of 
statistical physics in the last decades of the nineteenth 
century was to reconcile the time-asymmetric evolution 
of macroscopic systems with time-symmetric microscopic 
dynamics [l|. They came up with a probabilistic solution 



to this irreversibility paradox. Macroscopic states, they 
argued, are superpositions of an enormous amount of mi- 
croscopic states, the majority of them evolving in accor- 
dance with the second law of thermodynamics. The like- 
lihood that a macroscopic state violates the second law of 
thermodynamics is thus minute, typically exponentially 
small in the number of atoms it contains. Irreversibil- 
ity at the macroscopic level follows "by assuming a very 
improbable (i.e. with a very low entropy) initial state of 
the entire universe" [2j.|3J|. This mechanism works equally 
well in either quantum or classical systems. 

Simple mechanisms of irreversibility already exist at 
the microscopic level in chaotic (in particular mixing) 
classical systems with few degrees of freedom. As a 
matter of fact, mixing ensures that, after a sufficiently 
long evolution time, two initially well separated phase- 
space distributions will evenly fill phase-space cells of 
any given size. Since phase-space points can never be 
located with infinite precision, irreversibility sets in after 
mixing has occurred on a scale smaller than the phase- 
space resolution scale. This mechanism cannot be car- 
ried over to quantum systems, however, mostly because 
the Schrodinger time-evolution is unitary, in either real- 
or momentum-space. Microscopic quantum systems are 
generically stable under time-reversal, even when their 
classical counterpart is irreversible [J|. Peres instead sug- 
gested to investigate quantum irreversibility at the micro- 
scopic level through the fidelity 



M L (t) = KV'o \exp[iHt] exp[-iH t]\il>o)\' 



(1) 



with which a quantum state ipo can be reconstructed by 
inverting the dynamics after a time t with a perturbed 
Hamiltonian H = Hq + £ J£|. Because of its connec- 
tion with the gedanken time-reversal experiment pro- 
posed by Loschmidt in his argument against Boltzman's 
H-theorem nj|, Ml(£) has been dubbed the Loschmidt 



Echo by Jalabert and Pastawski |fj]. 

Echo experiments abound in nuclear magnetic reso- 
nance lit la, optics J2( , atomic [lfj , and condensed matter 
physics [ll|. Fundamentally, they are all based on the 
same principle of a sequence of electromagnetic pulses 
whose purpose it is to reverse the sign of the Hamilto- 
nian, H — > —Ho, by means of effective changes of co- 
ordinate axes J2J|. Imperfections in the pulse sequence 
result instead in H$ — > — Hq — S, and one therefore ex- 
pects the Loschmidt Echo to capture the physics of the 
experiments. This line of reasoning however neglects the 
fact that the time-reversal operation affects at best only 
part of the system, for instance because the system is 
composed of so many degrees of freedom, that the time 
arrow can be inverted only for a fraction of them. This is 
generically the case, as any system is coupled to an exter- 
nal, uncontrolled environment. To capture the physics of 
echo experiments one thus has to take into account that 
(i) the system decomposes into two interacting subsys- 
tems 1 and 2; (ii) the initial state of the controlled sub- 
system 1 is prepared, i.e. well defined, and its final state 
is measured and compared to the initial one; (iii) both 
the initial and final states of the uncontrolled subsystem 
2 are unknown; (iv) the Hamiltonian of system 1 is time- 
reversed with some tunable accuracy, however both the 
Hamiltonian of system 2 and the interaction between the 
two subsystems are uncontrolled. We therefore propose 
to investigate the physics of echo experiments by means 
of the following partial fidelity (we set h = 1) 



M B (f) = ((V'i|Tr 2 [ 



-i"K b t 



-i-Kit i-K t t & i-K^V\ 



|V>1>>,(2) 



where the forward and backward (partially time- 
reversed) Hamiltonians read 



Mb = -[Jfi + Ei]€ 



$H 2 + U f , 

h+h® [H 2 



(3a) 
S 2 ]+U b .(3b) 



The experiment starts with an initial density matrix 
Po = IV'iXV'il ® P2, which is propagated forward in time 
with %{. After a time t, we invert the dynamics of sys- 
tem 1. The imperfection in that time-reversal operation 



is modelled by Si, while £2 allows for system 2 to be 
affected by this operation (we will see below that trac- 
ing over the degrees of freedom of system 2 makes Mb (t) 
independent of either H 2 or £ 2 ). We leave open the pos- 
sibility that the interaction between the two systems is 
affected by the time-reversal operation, i.e. Iff may or 
may not be equal to lib. Because one has no control 
over system 2, the corresponding degrees of freedom are 
traced out. For the same reason, the outmost brackets in 
Eq. J2} indicate an average over p 2 ■ We name Mb (t) the 
Boltzmann echo to stress its connection to Boltzmann's 
counterargument to Loschmidt that time cannot be in- 
verted for all components of a system with many degrees 
of freedom. 

In this article, we present a semiclassical calculation 
of the Boltzmann echo for two classically chaotic subsys- 
tems along the lines of R.ef s . |(tL 1 1 Hi 1 1 4l | . and compare our 
results with those obtained from a Random Matrix The- 
ory (RMT) treatment of the problem. Our main result 
is that, in the regime of classically weak but quantum 
mechanically strong imperfection Si and coupling U^b, 
Mb(£) is the sum of two exponentials 

M B (i) =* cxp [- (r Sl + T f + T b ) t] + an cxp [-Ait] . (4) 

Here, ai is a weakly time-dependent prefactor, Ai is the 
classical Lyapunov exponent of system 1, and T^ and 
Tf.b are given by classical correlators for £1 and Uf,b 
respectively (see below). Equivalently, they can be re- 
garded as the golden rule width of the Lorentzian broad- 
ening of the levels of Hi induced by £1 and tl^b re- 
spectively |l2J. Together with the one- and two-particle 
level spacings 8i, 2 and bandwidths -Bi,2, they define 
the range of validity of the semiclassica l ap proach as 
5i < T El < Si, S 2 < r fjb < B 2 EBUU The sec- 
ond term on the right-hand side of Eq. <£jj exists exclu- 
sively for a classically meaningful initial state ipi such as 
a Gaussian wavepacket or a position state, but the first 
term is much more generic. It emerges from both a semi- 



classical or a RMT treatment and does not depend on 
the initial preparation tp\ of system 1. Other regimes of 
decay exist, which we here mention for the sake of com- 
pleteness. For quantum mechanically weak T^ <C <5i 
and Tfb <C 5 2 , one has a Gaussian decay, 



M B (*) = exp [- (£2/4 + U 2 /2 + Ug/2) V 



(5) 



in term of the typical squared matrix elements of £1 and 
Uf : b- Also, at short times a parabolic decay of Mb(£) 
prevails for any coupling strength. Finally, if system 1 is 
integrable, the decay of Mb(£) is power-law in time. 

The equivalence between Boltzmann and Loschmidt 
echoes is broken by Tf^, the decoherence rate of sys- 
tem 1 induced by the coupling to system 2 (or by U 2 h 
at weak interaction). Skillfull experimentalists can thus 
investigate decoherence in echo experiments with weak 
time-reversal imperfection £1 for which T^ 1 <C Tfb, and 
thus M B (i) ^ cxp[-(r f + r b )£] (or M B (t) ^ cxp[-(U| + 
U 2 ) t 2 /2] at weak interaction) as £1 is reduced. This 
might well be the explanation for the experimentally ob- 
served £i-independent decay of polarization echoes |"15l |. 

We now present our calculation. As starting point, 
we take chaotic one-particle Hamiltonians -ffi,2, and a 
smooth interaction potential U which depends only on 
the distance between the particles. We assume that 
it is characterized by a typical classical length scale, 
which in particular is larger than the de Broglie wave- 
length a of particle 1. For pedagogical reasons, we 
take narrow Gaussian wavepackets for the initial state 
of both particles, ifji(q) = (q\tpi) = (7rer 2 )~ di / 4 exp[ipi • 
(q — r.;) — |q — ri\ 2 /2cr 2 }. We note however that within 
our semiclassical approach, more general states can 
be taken for the uncontrolled system 2, such as ran- 
dom pure states p 2 = X^/3 a a a V\4>a){4>t3\, random mix- 
tures p 2 = ^2 a \a a \ 2 \(j) a )((j> a \ or thermal mixtures p 2 = 
J2 n exp [— f3E n ] \n)(n\. Arbitrary initial states for both 
subsystems can be considered within the RMT approach. 

From Eqs. l|2j and (jSJ we can rewrite Mb(£) as 



M B (t) = fdz 2 /j^dx^dq, Vi(qi)V'2(q2)V'l(q3)<q3,Z2|e- lMbt |xi,X2)(xi,X2|e- j:Hf *|qi,q2) 
J J i=i j=i 

We next introduce the semiclassical propagators (a = f, b labels forward or backward evolution; e") = — e 1 - 6 ' 
(xi,x 2 |e-^ t |qi,q 2 ) = £ 6^ cxp [i {e^>(qi, Xl ; t) + Sgfa, x 2 ; t) + §^ S2 (qi, x i; q 2 , x 2 ; t)} 



(6) 



(7) 



which are expressed as sums over pairs of classical trajectories, labeled Si (h) for particle i connecting q^ to X4 in the 
time t with dynamics determined by Hi or Hi + £j. Under our assumption of a classically weak coupling, classical 
trajectories are only determined by the one-particle Hamiltonians. Each pair of paths gives a contribution containing 
one-particle action integrals denoted by S Si (where we included the Maslov indices) and two-particle action integrals 
S«i,g 2 = / dr'Uf j b[qs 1 (T),q S2 (r)] accumulated along si and s 2 and the determinant 6 S] «, = C Sl C S2 of the stability 
matrix corresponding to the two-particle dynamics in the (di + d 2 )— dimensional space |lfil |. 



Our choice of initial Gaussian wave packets allows us to linearize the one-particle action integrals in q^ — r^ . We 

furthermore set Ssi,s 2 (qi, x i;q2, x 2;0 — §«i ,« 2 ( r i > x i i r 2, x 2! t), keeping in mind that r ± and r 2 , taken as arguments 
of the two-particle action integrals, have an uncertainty 0(a). We then perform six Gaussian integrations to get 

r 2 
M B (t) = (iira 2 ) 2 -^ I []dx i dy i dz 2 ]T A^^A^A^A^fC^ exp [i {$1 + $ 2 + $ 12 }] , (8) 



;=i 



paths 



where we wrote A Si = Ci i exp[— s 2 -(p Si — Pi) 2 ]- Paths with odd (even) indices correspond to system 1 (2). The 
semiclassical expression for M B (t) is obtained by enforcing a stationary phase condition on Eq. 0, i.e. keeping only 
terms which minimize the variation of the three action phases 



*i 



$2 = 



$12 



^(n.XiiO-Sf'o-i.xx;*) 
-Sj^n.yi.O+Si^ri.yi;*), (9a) 

^(r 2 ,x 2 ;i) + 4 b) ( x 2' z 2; i ) 

-^ ) (r2,y 2 ;t)-^ 4 b) ( y2 ,z 2 ;t), 

O(f) , o( b ) o(f) g(b) 

6 ai,«2 +%,* 2 - S i 3 ,«4 ~^3,U- 



(9b) 
(9c) 



The semiclassically relevant terms are identified by path 
contractions. The first stationary phase approximation 
over $i corresponds to contracting unperturbed paths 
with perturbed ones, Si ~ h and s 3 ~ l 3 . This pair- 
ing is allowed by our assumption of a classically weak 
Ei j]JD| • The phase $i is then given by the differ- 
ence of action integrals of the perturbation Si on paths 
si and s 3 , $i = 5S 8l (Ti,-x. 1 ;t) - 5S„ 3 (ri,yi,t), with 
SS Si = J Q drSi [q Si (r)] . Here, q Si (T) lies on Sj with 
q s ,(0) = ri and q Sl (t) = Xi, q S3 (t) = yi. A similar 
procedure for <& 2 requires s 2 ~ S4 and I2 — h, and thus 
x 2 ~ y 2 . These contractions lead to an exact cancellation 
$2 = 0, and one gets 

f 2 

M B (t) = (4 7 r ( j 2 ) H£l ^ /[]dx l dy J dz 2 ^(x 2 -y 2 ) 

i=l 

x £ |yi Sl | 2 |yL S2 | 2 K S3 | 2 |C i2 |e J [^-^3+*-]. ( io) 

Here, <5o-(x 2 — y 2 ) restricts the spatial integrations to |x 2 — 
y2 1 < a because of the finite resolution with which two 
paths can be equated. 

The semiclassical Boltzmann Echo ljlfl|l is dom- 
inated by two contributions. The first contri- 
bution is non diagonal in that all paths are 
uncorrelated. Applying the central limit the- 
orem one has (exp [i {5S Sl — SS S3 + $12}]) = 
exp [- (5SI) ((§ { ILr) <(S™ a ) 2 )] , where 

(^ 2 1 H/ t drdr'(I] 1 [q sl (r)]S] 1 [q sl (rO]}and<(§i f 1 ; b s ) 2 ) 2 ) 

= f*dr dr'(lX fjb [q sl (r),q 82 (r)] Uf )b [q 81 (r'),q» 2 (T / )]). 
In chaotic systems, correlators typically decay exponen- 
tially fast, thus (SS^) ~ T Sl t and ((S^) 2 ) ~ r f>b t. 



r 



Finally using the two sum rules 

(4™ 2 )^ /dx^|As,| 2 = L 

Si 

I dyj 8 a (y l - Xj) ^2 \ C U I = !' 



one obtains the nondiagonal contribution 

3v4 nd) (t) ~ cx P [-(r Sl +r f + r b )t]. 

The second contribution is diagonal, with s\ 1 
Xi ~ yi . From Eq. ifTTHl it reads 



(11a) 
(lib) 

(12) 
S3 and 



r(d) 



2d 1 +d 2 



M^(*) = (W)^ 






/ Y\ dx,dy,dz 2 S a (x.i - y,) 

i[AS. 1+ A8W , a +A 8 WJ ^ (13) 



X Z_^I s i K*-S2 1 |Cj 



where AS Sl = f Q drViSi[q sl (r)] • [q S3 (r) - q Sl (r)] and 

AS£ b) 2 = /odrViUf, b [q, I (r),q, !l (r)] • [q S3 (r) - q Sl (r)] . 
We perform a change of coordinates J dx.iJ2\C sl \ = 
/dpi, and use both the as vmp totics |C S1 | oc exp [—Ait] 
valid for chaotic systems |ljj and the sum rules of 
Eqs. (ffTTfl to get 



M^ d) (i) ~ ai exp [-Ait] 



(14) 



Here, ot\ is only algebraically time-dependent with a\(t = 
0) = 0(1). Together, diagonal ((141) and nondiagonal 
ltl2)l contributions sum up to our main result, Eq. Q. 
We finally note that the long-time saturation at the in- 
verse Hilbert space size of system 1, Mb (00) = N^ , 
is obtained from Eq. JHJ with the contractions si ~ S3, 
s 2 ~ S4, l\ ~ I3 and I2 — h- 

Analyzing Eq. J3J, we first note that M B (t) depends 
neither on i? 2 nor on S 2 . This is so because one traces 
over the uncontrolled degrees of freedom. We stress that 
this holds even for classically strong E 2 . Most impor- 
tantly, besides strong similarities with the Loschmidt 
Echo, such as competing golden rule and Lyapunov de- 
cays [y, |l2|, the Boltzmann Echo can exhibit a £j.~ 
independent decay given by the decoherence rates Ff jb 




Figure 1: Main plot: Boltzmann echo for N = 1024, 
A"i = K 2 = 10.09, and <n = 0.0018 (P Sl ~ 0.09). Data 
have been calculated from 50 different initial states. The full 
lines correspond to e = 0, 0.0018 and 0.0037 (from right 
to left) and the dashed lines give the predicted exponen- 
tial decay of Eq. Q), with T u = 1.210 4 e 2 ,r El = 2.610V?, 
A = 1.6 > ru,Tsi (dashed lines have been slightly shifted 
for clarity). The dotted line gives the saturation N~ . Inset : 
M B for e = 0.0037, and <n = 0.0003 (circles; T Sl ~ 2. 10~ 3 ), 
ax = 0.0006 (squares; T El ~ 9. 10 3 ), and 0.0009 (diamonds; 
F El — 0.02). The dashed line indicates the theoretical pre- 
diction Mb(£) = exp[-0.3t]. 



in the limit T^ 1 <C T ih . Extending our analysis to the 
regime T^ 1 <C 5\, Tfb <C $2 by means of quantum pertur- 
bation theory, we find a gaussian decay of Mb (t) , Eq. (§J . 
It is thus possible to reach either a Gaussian or an expo- 
nential, Ei-independent decay, depending on the balance 
between the accuracy Si with which the time-reversal op- 
eration is performed and the coupling between controlled 
and uncontrolled degrees of freedom. This might explain 
the experimentally observed saturation of the polariza- 
tion echo as £1 is reduced (jjji though a more precise 
analysis of these experiments in the light of the results 
presented here is necessary. 

We numerically illustrate our findings. We consider 
two coupled kicked rotators with Hamiltonian 

Hi = tfft + Kicmfa) ^2d(t-nT), (15a) 

n 

U = e sin(xi - x 2 - 0.33) ^ S(t - nT). (15b) 

n 

We concentrate on the regime Ki > 7, for which the 
dynamics is fully chaotic with Lyapunov exponent A, ~ 
ln[ifi/2]. The time- reversed one-particle Hamiltonians 
are obtained through Ki — ► Ki+<n. We here restrict our- 
selves to the case IX = Uf = lit,. Both rotators are quan- 
tized on the torus with discrete momenta p n = 2irn/N, 
n = 1,2, ...N. The one- and two-particle bandwidths 
and level spacings are given by B± = 2n, Si = 2tt/N and 
B2 = 47T, 82 = 4tt/N 2 . For more details on the numerical 
procedure, we refer the reader to Ref. |lS| . 

We first checked that Mb (t) is independent of K 2 (as 



long as system 2 remains chaotic) and a 2 , and therefore 
set K 2 = K\, (T2 = 0. The main panel in Fig. ^ shows 
that for Bi > T Sl > S 1: B 2 > T u > S 2 , Eq. (2J is sat- 
isfied. Additionally, the inset of Fig. Q] illustrates that 
when r^j -C 2Fu, the observed decay is only sensitive to 
U , and one effectively obtains a £i-independent decay. 
Further unshown data confirm the existence of the Lya- 
punov decay [second term in Eq. <Q}]. All our numerical 
results thus confirm the validity of Eq. (QJ . 

In conclusion we propose to analyze echo experiments 
in the light of the Boltzmann echo of Eq. (2J and 10 ■ 
Our semiclassical and RMT analysis showed that the de- 
cay of Mb(£) saturates at a finite value even when the 
time-reversal operation is performed with infinite accu- 
racy. Further work should attem pt t o connect these re- 
sults with echo experiments |sLldllfl llllll.'i| . 
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